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ABSTRACT 

Context. A debated topic in star formation theory is the role of magnetic fields during the protostellar phase of high-mass stars. It 
is still unclear how magnetic fields influence the formation and dynamics of massive disks and outflows. Most current information 
on magnetic fields close to high-mass protostars comes from polarized maser emissions, which allows us to investigate the magnetic 
field on small scales by using very long-baseline interferometry. 

Aims. The massive star-forming region W75N contains three radio continuum sources (VLA 1, VLA2, and VLA3), at three different 
evolutionary stages, and associated masers, while a large-scale molecular bipolar outflow is also present. Very recently, polarization 
observations of the 6.7 GHz methanol masers at milliarsecond resolution have been able to probe the strength and structure of the 
magnetic field over more than 2000 AU around VLA 1. The magnetic field is parallel to the outflow, suggesting that VLA 1 is its 
powering source. The observations of H 2 masers at 22 GHz can give more information about the gas dynamics and the magnetic 
fields around VLA 1 and VLA 2. 

Methods. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz 
water masers in the star-forming region W75N. 

Results. We detected 124 water masers, 36 around VLA 1 and 88 around VLA 2 of W75N, which indicate two different physical 
environments around the two sources, where VLA 1 is in a more evolved state. The linear polarization of the masers confirms the 
tightly ordered magnetic field around VLA 1 , which is aligned with the large-scale molecular outflow, and also reveals an ordered 
magnetic field around VLA 2, which is not parallel to the outflow. The Zeeman-splitting measured on 20 of the masers indicates 
strong magnetic fields around both sources (the averaged values are |BvlaiI ~ 700 mG and |Bvla2| ~ 1700 mG). The high values of 
the magnetic field strengths, which come from the shock compression of the gas, are consistent with the methanol and OH magnetic 
field strengths. Moreover, by studying the maser properties we were also able to determine that the water masers are pumped in 
C-shocks in both sources. 

Key words. Stars: formation - masers: water - polarization - magnetic fields - ISM: individual: W75N 

1 . Introduction They obtained their results by using methanol masers as probes 

of magnetic fields. In fact, the bright and narrow spectral line 
During the formation of low-mass stars, the magnetic field is emission of masers is ideal for detecting, with polarimetric in- 
thought to slow the collapse, to transfer the angular momentum, terferometry, both Zeeman-splitting and the orientation of the 
and to power the outflow, but its role during the formation of ma g ne tic field on scales from arcsec to milliarcsec (mas). So far, 
high- mass stars i s still under debate (e.g., McKee & Ostriker me main maser species used for this purpose are H 2 0, OH, and 
w , |2007J Girart et al. |2009J. Several questions remain unanswered C H 3 OH masers. Since their masing conditions are different they 
because massive star-forming regions are rare and distant and of- are located in distinct zones f massive star-forming regions, 
ten consist of a large number of protostars, consequently it is dif- xhe Hz0 masers are det ected in the denser zones with hydro- 
ficult to observe magnetic fields during their protostellar phase. gen num ber densities n H2 between approximately 10 8 and 10 10 
Recently some progress has been made. Vlemmings et al. (2010) cm -3 (Elitzur et al. [19591 . showing a velocity width of about 
and Surcis et al. (|2009[ hereafter S09) have shown that the mag- 1 s -i and a high brightness temperature T h > 10 9 K (e.g., 
netic field is orthogonal to large rotating disks and parallel to Reid & Moran 1981 , xhe first very long baseline interferome- 
molecular bipolar outflows even during the high-mass protostel- try (V LBI) observations of the linear and circular polarization of 
lar phase and not only during the fo rmation of low-mass stars water maser emissio ns were made by Leppanen et al. JT9981 and 
(e.g., Matsumoto & Tomisaka|200l McKee & Ostaker|2007J. Sarma et al. (I20011 . respectively. Afterwards more H 2 maser 

polarization observations were carried out by various authors, 

* Member of the International Max Planck Research School (IMPRS) which confirmed the importance of this kind of observation in 
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(e.g., Imai et al. 120031 Vlemmings et al. 12006b . Here we present 
Very Long Baseline Array (VLBA) linear and circular polariza- 
tion observations of water masers in the high-mass star-forming 
region W75N. 

W75N is an active high-mass star-forming region in the 
molec ular complex DR21-W75 (Dickel et al. 119781 Persi et al. 
2006) at a distance of 2kpc. At the resolution of ~ 175 three 
Hn regions were identified (Haschick et al. |19811 >: W75N(A), 
W75N(B), and W75N(C). At 075 resolution, Hunter et al. 
(1994) detected three subregions in W75N(B): Ba, B b and Be. 
At higher angular resolution (~ 071), Torrelles et al. (1997]) re- 
solved the subregions Ba and Bb further (renaming them VLA 1 
and VLA 3), and imaged a third weaker and more compact H n 
region between them, which they named VLA 2. A large-scale 
high-velocity outflow, with an extension greater than 3 pc and 
a total molecular mass greater than 255 M , was also detected 
from W 75N (B) (e.g., Shepherd et al. I2001T ). Shepherd et al. 
(2003) proposed a multi-outflow scenario where VLA 2 may 
drive the large-scale outflow, and VLA 1 and VLA 3 are the 
centers of two other small flows. So far, it has been impossi- 
ble to determine the main powering source of the 3 pc outflow. 
Several authors have suggested VLA 1 as the powering source 
(e.g., Torrelles et al. 1997 S09 and references therein). 

Several maser species (H2O , OH, and CH3OH) have been 
detected in W75N (e.g., Torrelles et al. [19971 Baart et al[1986 
S09), in particular around the two H 11 regions VLA 1 and VLA 2. 
The H2O maser emissions associated with VLA 1 are located 
along a linear structure parallel to its radio jet with a proper 
motion of about 2 mas/yr, while those associated with VLA 2 
show a circular distribution, which is expanding with a veloc- 
ity of about 5 mas/yr (Torrelles et al., 120031 hereafter T03). 
Only one H2O maser is associated with VLA 3 (T03). Methanol 
masers were detected in two groups (A and B/C) located north- 
west and southeast of VLA 1, respectively. Group A is along a 
linear structure and group B/C show an arc-like structure (S09). 
No CH3OH masers are associated with VLA 2 (e.g., Minier et al. 
2000 S09). However, VLA 2 is the place where the most inten- 
sive OH flare took place (Alakov et al. 120051 Slysh et al. EOTOl ). 
Other OH maser emission sites are situated on a ring structure 
around all three radio sources (Hutawarakorn et al. 2002). Based 
on the different activity in H2O and OH masers, Torrelles et al. 
(119971 ) suggest that these three sources are at different evolu- 
tionary stages, in particular VLA 1 is the "oldest" source of this 
region and VLA 3 the "youngest" one (driving a radio Herbig- 
Haro object, Carrasco-Gonzalez et al. 2010). In this evolutionary 
scenario, VLA 2 is at an intermediate stage. 

In order to determine and investigate the magnetic field of 
W75N(B), polarization observations of OH and CH3OH maser 
emissions have been made (e.g., Hutawarakorn et al. [2002; 
Fish & Reid 120071 S09). The magnetic field strength obtained 
from the OH maser polarization observations was ~1 mG (e.g., 
Hutawarakorn et al. 120021 Slysh et al. 120021) . Dur ing the OH 
maser flare near VLA 2, Slysh & Migenes (2006) detected a 
strong magnetic field of more than 40 mG in several maser spots, 
which increased in the next observations up to 70 mG (Slysh et 
al. 120101 ). S09 measured a magnetic field strength of 50 mG 
around VLA 1 by studying the circular polarized emissions of 
methanol masers. Hutawarakorn & Cohen (1996) and S09 in- 
vestigated the linear polarization of OH and CH3OH masers, re- 
spectively, and both suggested that the magnetic field is oriented 
along the outflow. Although the OH masers at 1.6 and 1.7 GHz 
are very susceptible to both internal and external Faraday rota- 
tion, at 6.7 GHz the Faraday rotation is only about 5° indicating 
that the finding of the large-scale magnetic field aligned with the 



outflow is robust. 

Since the water masers, unlike the CH3OH and OH masers, 
are found in denser zones of the star-forming regions, it is crucial 
to study their polarization emission on a small scale. In fact, be- 
cause they arise close to both radio sources they can give us new 
information on the role of the magnetic field during the massive 
star formation, in particular at two different evolutionary stages. 

2. Observations and data reduction 

We observed the massive star-forming region W75N(B) in the 
6i6-523 transition of H2O (rest frequency: 22.23508 GHz) with 
the NRAC0 VLBA on November 21 st 2005. The observations 
were made in full polarization spectral mode using 4 overlapped 
baseband filters of 1 MHz in order to cover a total velocity range 
of « 50 kms _I . Two correlations were performed. One with 
128 channels in order to generate all 4 polarization combina- 
tions (RR, LL, RL, LR) with a spectral resolution of 7.8 kHz 
(0.1 kms 1 ). The other one with high spectral resolution (512 
channels; 1.96 kHz=0.027 kms 1 ), which only contains the cir- 
cular polarization combinations (LL, RR), to be able to detect 
Zeeman splitting of the H2O maser across the entire velocity 
range. Including the overheads, the total observation time was 
8h. 

The data were reduced using the Astronomical Image 
Processing Software package (AIPS) following the method of 
Kemball et al. ([T995). The bandpass, the delay, the phase 
and the polarization calibration were performed on the cal- 
ibrator J2202+4216, which has been used successfully by 
Vlemmings et al. (2006 ) for H2O maser polarization observa- 
tions of CepheusA. The fringe-fitting and the self-calibration 
were performed on one of the brightest maser features VLA 1.14 
(Table [TJ. All calibration steps were initially performed on the 
dataset with modest spectral resolution after which the solu- 
tions, with the exception of the bandpass solutions that were 
obtained separately, were copied and applied to the high spec- 
tral resolution dataset. Stokes / (rms - 7.3 mJybeanT 1 ), Q 
(rms = 6.5 mJybearrT 1 ) and U (rms = 6.5 mJybeanT 1 ) data 
cubes were created using the AIPS task IMAGR (beam-size 
2.0 mas x 0.7 mas) from the modest spectral resolution dataset, 
while the / and V (rms =10 mJybeairT 1 ) cubes were imaged 
from the high spectral resolution dataset and for the same fields. 
The Q and U cubes were combined to produce cubes of polarized 
intensity and polarization angle. Since these observations were 
obtained between two VLA polarization calibration observa- 
tionfl during which the linear polarization angle of J2202+4216 
was constant at 13?5, we were able to estimate the polarization 
angles with a systemic error of no more than ~2°. The / and 
V cubes at high spectral resolution were used to determine the 
Zeeman splitting. 

3. Analysis 

3.1. Maser identification 

To identify the water maser features in W75N(B) region, we 
used an identification process divided into three steps: 1) we 
used a program called "maser finder" (Curiel, personal com- 
munication), which is able to search for maser spots, velocity 

1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

2 http://www.aoc.nrao.edu/astro/calib/polar 
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Fig. 1. Results of the full radiative transfer^ 2 -model fits, for the maser that shows the highest linear polarization fraction (VLA1.19) 
and for the maser that shows the highest Py around VLA2 (VLA2.22). The fits yield the emerging maser brightness temperature 
7bAQ and the intrinsic maser linewidth AV\. Contours indicate the significance intervals Ay 2 =0.25, 0.5, 1, 2, 3, 7, with the thick 
solid contours indicating lcr and 3<x areas. 






Fig. 2. Total intensity and circular polarization spectrum for the three maser features VLA1 .07, VLA1 . 16, and VLA2. 16. The thick 
red lines are the best-fit models of I and V emission obtained using the radiative transfer method (see Sect. 13.21 ). All maser features 
were centered to zero velocity. 



channel by velocity channel, with a signal-to-noise ratio greater 
than a given value (in our case, an 8-sigma limit was used); 2) 
we fitted the identified maser spots using the AIPS task IMFIT; 
and 3) we identified a maser feature, when three or more maser 
spots coincided spatially (within a box of 2 by 2 pixels) and each 
of them appeared in consecutive velocity channels. Table[T]con- 
tains the positions, velocities and peak flux densities of all the 
maser features found in the region. 

In the first part of the process and after evaluating the global 
noise level of the channel map, the code searches for maser spots 
in the velocity map, excluding the edges of the channel map. 
When strong maser spots are found (for instance, with a peak 
flux density higher than 2 Jy beam 1 ), the local rms around these 
strong maser spots (within a predefined box, for instance, 400 
by 400 pixels) is estimated, as well as the rms outside the boxes. 
The code identifies a maser spot when the SNR of the candi- 
date (using the local rms) is greater than a predefined value. In 
the case of W75N, we have adopted a lower limit of 8-sigma. 
The code produces a table with the spatial location and intensity 
of all the identified maser spots. This procedure is repeated for 
each velocity channel included in the input file that is used by the 
code. The code also produces a script for each velocity channel 
that can be run inside AIPS in order to produce Gaussian fits of 



all the identified maser spots. After running all the AIPS scripts, 
a table with all the important parameters (e.g., spatial resolution, 
velocity, intensity) of each maser spot is produced. We then use 
the output table to search for maser features that fulfill the crite- 
ria described above. 

3.2. Polarization fitting 

From the maser theory we know that the fractional linear po- 
larization Pi of the H2O maser emission depends on the degree 
of its saturation and the angle (9) between the maser propaga- 
tion direction and the magnetic field (e.g., Goldreich et al. 119731 
Nedoluha & Watson [1992| l. Even the relation between the mea- 
sured polarization angle x an d the magnetic field angle on the 
sky (<Pb) depends on 8, with the linear polarization vector per- 
pendicular to the field for 6 > cl ; t ~ 55° (Goldreich et al. 
119731 ). where cr j t is the so-called Van Vleck angle. It is there- 
fore important to evaluate 6 for every H2O maser that shows 
linear polarization emission. In order to obtain the best values 
for 0, we needed to determine the emerging brightness tempera- 
ture (7bAQ) and the intrinsic thermal linewidth (AVj), which is 
the full width half-maximum (FWHM) of the Maxwellian distri- 
bution of particle velocities, of the masers with high precision. 
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Table 1. All 22 GHz water maser features in W75N for which linear and/or circular polarization was detected. 
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cl The water masers associated with VLA 1 are divided in three groups according to their positions (see Sect. l4.ll . 

h The best-fitting results obtained by using a model based on the radiative transfer theory of water masers (Vlemmings et al. 2006 2010 1 for Y + Y v = 1. The errors were determined by 
analysing the full probability distribution function. For T ~ 400 K (!"„ = 2) and T ~ 2500 K (r„ = 13) r b A£2 has to be adjusted by adding +0.48 and +1.15, respectively (Nedoluha & 
Watson fT99"2l Anderson & Watson f!993T 

c The angle between the magnetic field and the maser propagation direction is determined by using the observed Pi and the fitted emerging brightness temperature.The errors were determined 
by analysing the full probability distribution function. 

''Because of the degree of the saturation of these water masers 7bA£l is underestimated. AV; and are overestimated. 

The constraint fit did not allow us to evaluate the error bars properly. 

^In the fitting model we include the values for 7\, AH and A V\ of the closest features. 



This could be done by fitting the water maser emissions with 
the full radiative transfer method described by Vlemmings et al. 
(20061 120101) . who successfully used it for water and methanol 
masers in Cepheus A. These radiative transfer methods are based 
on the models for water masers of Nedoluha & Watson ( 1992 ), 
for which the shapes of the total intensity, linear polarization, 
and circular polarization spectra depend on 7t,Af2 and AVj. They 
also found that the emerging brightness temperature scaled lin- 



early with (T + T v ), which are the maser decay rate T and cross- 
relaxation rate F v . Since for the 22 GHz H2O masers T < Is -1 
and r y depends on the gas temperature, we include in our fit a 
value (F + F v ) = Is -1 , so that we could adjust the fitted 7b AO 
values by simply scaling according to (r + T v ). Note that AV; and 
9 do not need to be adjusted. 

We model the observed linear polarized and total intensity 
maser spectra by gridding the intrinsic thermal linewidth AV,- 
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between 0.4 and 3.5 kms -1 , in steps of 0.025 kms -1 , using a 
least square fitting routine (x 2 -model) with an upper limit of 
the brightness temperature T\,AQ. - 10 11 K sr. In Fig. [T]two 
examples are shown of the results obtained by fitting the wa- 
ter maser emissions with the full radiative transfer method code. 
From these first fits, we can then obtain the angles 9 by consid- 
ering the relation between Pi and 9 (see Vlemmings et al. 120101 
for more details). The best values for T^AQ. and AVj are then in- 
cluded in the full radiative transfer code to produce the / and V 
models that were used for fitting the total intensity and circular 
polarized spectra of the H2O masers extracted from the high- 
resolution / and V cubes. The magnetic field strength along the 
line of sight is evaluated by using the equation 



B11 = B cos9 



Py Av L 
2-A 



(1) 



F-F' 



where Avi is the FWHM of the total intensity spectrum, the 
Ap-F coefficient, which depends on T^ACl. The circular polar- 
ization fraction (Py) are obtained from the / and V models. 

Another important aspect that must be considered in the 
analysis of the water masers is the degree of their saturation, 
which influences Pi. We are able to estimate it by considering 
the ratio between maser rate of stimulated emission (R) and the 
sum of maser decay rate and the cross-relaxation rate (F + T v ). 

The stimulated emission rate is given by 



R 



Ak B T b A£l 
Anhv 



(2) 



where A = 2 x 10~ 9 s _1 is the 22-GHz H2O maser spontaneous 
rate (Goldreich & Keeley [T972l l. ks and h are the Boltzmann and 
Planck constants, respectively, v the maser frequency and T^ACl 
the emerging brightness temperature adjusted at the right F + T v 
value. The masers are unsaturated when R/(T + F v ) < 1, in the 
onset of saturation when R / (F+r v ) > 1 , and fully saturated when 
R/(T + T v ) =3 100 (Vlemmings et al. l20061 >. Moreover, when the 
saturation sets in, the maser lines start to broaden again, and this 
means that the observed linewidth (Avl) becomes as large as AVj 
or larger depending on the degree of the saturation. 

Although the full radiative transfer method is based on a 
model for general water masers, it should be noted that, when 
this model is applied to saturated masers, it is impossible to 
properly disentangle the values of T\,AD. and AV;. In this case 
the method only provides a lower limit for r D AQ and an up- 
per limit for AVj. For high maser brightness temperatures, i.e. 
greater than 10 9 crrT 3 K sr (Nedoluha & Watson 1992), the cosf? 
dependence of Eq. Q] breaks down thereby introducing a more 
complex dependence on 9. As a result, the angle 9 obtained from 
the fit is overestimated and, in particular, it could also reach val- 
ues of 90° or greater. Consequently, the values of T^ACl, AVi and 
9 obtained for saturated masers could not be taken into account 
in the analysis of the region. 



4. Results 

4.1. Maser distribution 

In Fig. [3] we show the distribution of the 124 22-GHz water 
maser features detected with the VLBA. No water maser emis- 
sion with a peak flux density less than 80 mJybeairT 1 is de- 
tected. All H2O masers are detected around the radio sources 
VLA 1 (29%, 36/124) and VLA2 (71%, 88/124). 

The water masers associated with VLA 2, which have local 
standard of rest velocities (Vlsr) between -7.7 and 13.7 kms -1 , 



are distributed elliptically around its 1.3 cm continuum peak. In 
Fig. 3] we report two ellipses obtained by fitting the positions 
of the water masers detected by T03 (called here ellipse 1) and 
in the present paper (ellipse 2), their parameters are reported in 
Table [2] The major axis of ellipse 2 is 64 mas larger than that 
of ellipse 1, the minor axes are also in the same ratio. This in- 
dicates an expansion velocity of ~4.8 mas/yr in all directions, 
corresponding to ~ 46 kms -1 (at 2 kpc). This value is con- 
sistent with the proper motion of the water masers as reported 
by T03. Since the expansion velocity is symmetric in all direc- 
tions, we aligned the centers of the two ellipses with the posi- 
tion of the 1.3 cm continuum peak obtained from a Gaussian 
fit. Consequently we were able to overlay all water masers of 
the region on the continuum VLA map with an accuracy of ap- 
proximately 10 mas. In this way we could also compare the po- 
sitions of the masers around VLA 1 at the two different epochs. 
If we look at the northern masers around VLA 2, we see that the 
masers detected by T03 are better aligned with the ellipse 1 than 
the masers detected by us with ellipse 2 (Fig. |4] right panel), in- 
dicating that those masers are also moving northeastward. This 
movement might be due to the formation of a jet (see Fig . [3j . 

The water masers associated with VLA 1 can be divided into 
three groups, one composed of the northern masers (VLA 1A), 
one of the masers closest to the central protostar (VLA IB, see 



Table 2. Number, epoch, semi-major axis, semi-minor axis, and 
inclination from the elliptical fit of the water maser associated 
with VLA 2 in the two different epochs. 



Ellipse 


epoch 


a 

(mas) 


b 

(mas) 


8 

(°) 


ref. 


1 

2 


2-Apr.-1999 
21-Nov.-2005 


71.0 
103.1 


61.1 

92.9 


123.5 
120.5 


T03 
this work 



Fig. [5] right panel) and the third one of the southern masers 
(VLA 1C). The east-northern masers of the VLA 1A group did 
not show emission during the previous observation of epoch 
1999 (T03). All VLA 1A masers (Vlsr - 19-25 kms" 1 ) 
are well-aligned with the direction of the thermal radio jet de- 
tected by Torrelles et al. d 19971 1 and the methanol masers de- 
tected by S09. The VLA IB and VLA 1C groups show a dis- 
tribution quite similar to that shown on April 2 nd 1999 (T03). 
For these two groups T03 reported a proper motion of about 
~2 mas/yr (~ 19 kms -1 ). Although we have detected several 
water masers close to those detected by T03, it is very difficult 
to estimate the proper motion between the two epochs. In fact, 
as is shown on the right panel of Fig. [5] it is difficult to iden- 
tify a maser that arises in both epochs (1999 and 2005). The 
LSR radial velocities of VLA 1 are between 9 and 25 kms -1 , 
which are about 3-4 times higher than those of methanol masers 
(S09) and greater than the systemic velocity of the whole region 
(10 kms 1 . Shepherd et al. 120051 . 

4.2. Linear and circular polarization 

Linear polarization was detected in 50 masers (Table[T]i. In par- 
ticular, about 70% of the H2O masers associated with VLA 1 
show linear polarization emission. The corresponding linear 
polarization vectors are shown in Figs. |4] and [5] (left panel). 
The highest fractional linear polarization was detected in wa- 
ter masers related to VLA 1 (Pi ~ 26%). Eight masers associ- 
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Fig. 3. Positions of water and methanol masers superimposed on the 1.3 cm continuum contour map of the VLA1 thermal jet and 
VLA 2 (Torrelles et al. 1997). Contours are -3, -2, 2, 3, 4, 5, 6, 7, 8 x 0.16 mJybeairT 1 . Blue triangles indicate the positions of 
the water maser features detected with the VLBA (this paper). Red circles indicate the position of the methanol masers features 
detected by S09. The two arrows indicate the direction of the large-scale molecular bipolar outflow (PA=66°). 



ated with VLA 1 have high linear polarizations, Pi between 10% 
and 26%, so that the median value of linear polarizations for the 
masers in VLA 1 is ~7%. In contrast, the maximum Pi for the 
masers associated with VLA 2 is 6%, and the median value is 
much lower at 1.5%. Circular polarization was detected in 20 
masers, equally distributed between VLA 1 and VLA 2. 

The rms weighted linear polarization angles of the water 
masers around VLA 1 and VLA 2 are (xi) ~ -67°+40° and 
(X2) ~ -72 ±32°, respectively. Since the water masers of 
group VLA 1A show ordered linear polarization vectors, it is 
worthwhile evaluating its weighted polarization angles, which 
is Ctia) = -20° ±9°. 

We were able to fit 42 masers with the method described 
in § 13.21 and the results are given in columns 10 and 11 of 
Table [TJ and two examples (i.e., VLA1.19, VLA2.22) are re- 
ported in Fig.Q] Although the emerging brightness temperatures 
for VLA 1 and VLA 2 are both in the range 10 8 Ksr < T b AQ. < 
10 10 K sr, their weighted intrinsic maser linewidths are very dif- 
ferent. They are (AV^vlai = LOtg® kms" 1 and <AVi) V LA2 = 
respectively. 

Finally, considering the emerging brightness temperature 
and the observed Pi we determined 9 for those masers which 
we were able to fit for AVj and 7bA£2. The results are given 
in column 14 of Table Q] with the errors determined by ana- 



2.5t°ol km s" 1 



lyzing the full probability distribution functions. Eight of the 
H2O masers around the radio source VLA 1 (from VLA 1.1 2 to 
VLA1.19, which are the closest features to the protostar) show 
values for 9 equal to 90°. This suggests that, from the theory of 
the Zeeman effect, we are observing perpendicular to the mag- 
netic field. Consequently, we should not be able to detect circu- 
lar polarization from these masers. However, since we detected 
circular polarization in three of them (VLA1.16, VLA1.17, and 
VLA1 . 19) the measurements of 9 must be affected by the degree 
of their saturation, as discussed in § 13.21 By excluding these 8 
masers, we evaluated <0> V lai = 83°+7° and (0) VL A2 = 85°+6°, 
which imply that the magnetic fields are close to the plane of the 
sky. 

In column 13 of TableQ]we report the magnetic field strength 
along the line of sight obtained by fitting the high-resolution 
spectra of / and V (see §E3. The masers VLA1.21, VLA2.01, 
and VLA2.02 do not show any linear polarization emission, but 
they do show circular polarization emission. Consequently we 
could not get TbAQ and AVj as done for the other masers that 
showed linear polarization emission, and to measure the mag- 
netic field strength, we used the T\,A£l and AV; values of the 
closest masers. In Table Q] the circular polarization fraction is 
also reported (column 12). 
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Fig. 4. Left panel: a close-up view of the H2O maser features around the radio source VLA2. Right panel: a zoom-in view of 
the boxed region of the left panel. The octagonal symbols are the identified maser features in present work scaled logarithmically 
according to their peak flux density. The maser LSR radial velocity is indicated by color. A 5 JybeairT 1 symbol is plotted for 
illustration in both panel. The linear polarization vectors, scaled logarithmically according to polarization fraction Pi (in Table[TJ, 
are overplotted. Two ellipses are also drawn in both panels. They are the results of the best fit of the water masers (crosses) detected 
by T03 (dotted ellipse, ellipse 1; epoch 1999) and of those detected in present work (solid ellipse, ellipse 2; epoch 2005). Their 
parameters are listed in Table|4] The synthesized beam is 2.0 mas x 0.7 mas. 
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Fig. 5. Left panel: a close-up view of the H2O maser features (octagon) around the radio source VLA 1 . Right panel: a zoom-in view 
of the water masers group VLA IB, here are also reported the water masers (square) detected by T03. The octagonal and the square 
symbols are scaled logarithmically according to their peak flux density. The maser LSR radial velocity is indicated by color. The 
systemic LSR radial velocity of W75N massive star-forming region is 10 kms -1 (Shepherd et al. I2003I I. Symbols for 3 Jybeam -1 
and for 1 Jy beam 1 are shown in the lower left corner of the panels. The linear polarization vectors, scaled logarithmically according 
to polarization fraction Pi (in Table [TJ, are overplotted on the left panel. The synthesized beam is 2.0 mas x 0.7 mas. 



5. Discussion 



The water masers detected using the VLBA are related to two 
different radio continuum sources of W75N(B), VLA 1 and 
VLA 2, which are thought to be in two different evolutionary 
stages. In the next sections, we discuss their nature as two inde- 



pendent groups, their polarization and finally the type of shocks 
that give rise to them. 
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5.1. H 2 maser properties 

As suggested by the difference between their intrinsic thermal 
linewidths, which is 1.5 kms -1 , the masers around VLA 1 and 
VLA2 appear to arise under two different physical conditions. 
Using 

AV; * 0.5 X (77100)5, (3) 

we estimate the gas temperatures in the H2O masing regions 
around VLA 1 and VLA 2, and find them to be TVla 1 ~ 
400 K and 7\/la2 ~ 2500 K, respectively. The intrinsic ther- 
mal linewidths of the VLA 2 masers are quite similar to each 
other, while the three groups of masers around VLA 1 show 
a difference in their AVj. In particular, we find that the maser 
groups VLA 1A and C have similar intrinsic thermal linewidths 
((AV;) = 1.0 kms -1 ) while the group VLA IB has lower values 
((AV;) = 0.7 kms -1 ). This can come from the degree of their 
saturation that we can evaluate by considering Eq.|2] 

Before using Eq. [2] we have to adjust the emerging bright- 
ness temperature values as explained in § 13.21 At a temperature 
of -400 K, for which T + r v = 3 s 1 (Nedoluha & Watson [T992l , 
the weighted emerging brightness temperature of VLA IB is 
(TbAO) > 10 11 Ksr, which is close to the limit of full satura- 
tion, and for VLA 1A and C is <r b AO> = 10 1() K sr. From Eq.|2] 
the masers of VLA 1A and VLA 1C are thus completely unsat- 
urated (R/(T + r v ) = 0.5), while VLA IB masers are saturated 
(R/{r + F v ) > 5), even if the real degree of their saturation is un- 
known because of underestimating TbAO. The average observed 
linewidth (Av'l) of VLA IB group, which is about 0.5 kms -1 , 
gives us another indication of the saturation state of their H2O 
masers. As reported in § 13.21 the intrinsic thermal linewidths for 
the saturated masers are overestimated so Avl might be larger or 
at least equal to AV; implying that the maser lines are rebroad- 
ened as expected in the saturated maser lines. As the VLA IB 
masers are saturated, they have not been taken into account in 
our conclusions. For the other masers associated with VLA 1 , we 
observed linewidths well below their intrinsic thermal linewidth, 
and this confirms that these masers are unsaturated. The masers 
associated with VLA 2 are unsaturated; in fact, for T ~ 2500 K 
(r + T v = 14s~\ Anderson & Watson [19931 1. we get ratios of 
fl/(r + F v ) <0.5. 

From the measurements of the maser flux densities (S(v)) 
and feature angular sizes (£) we can estimate the brightness tem- 
perature by the equation reported in S09, 

T b _ S(y) ( 2 2 r 1 

where £h 2 o = 1-24 x 10 9 mas 2 Jy~' K is a constant factor which 
includes all constant values, such as the Boltzmann constant, 
the wavelength, and the proportionality factor obtained for a 
Gaussian shape by Burns et al. ( 1 19791 1. Comparing 7b with the 
emerging brightness temperatures obtained from the model, we 
can easily estimate AO. The Gaussian fit of the masers give a 
size of 0.4 mas for the water masers related to VLA1, which 
indicates that they are unresolved, while a size of 0.5 for the 
masers associated with VLA 2, which indicates that they are 
marginally resolved. The brightness temperatures of the bright- 
est masers in both groups are 7i,,vlai.08 > 1-23 x 10 12 K and 
7b,VLA2.24 ~ 10 12 K, for which we find AQvlai.08 < 10~ 2 sr 
and AQ VL A2.24 ~ 10~ 2 sr for T = 400 K and 2500 K, respec- 
tively. Considering all the features in Table [1] the maser beam- 
ing of the two groups are AOvlai ^ 10~ 2 sr (T + T v = 2s _1 ) 



and A£Xvla2 ~ 10~ 2 sr (r + F v = 14 s~'). In a tubular geom- 
etry AQ w (d/l) 2 , where d and / are the transverse size and 
length of the tube respectively. Assuming d approximately the 
size of the maser features, the maser lengths are all in the range 
10 12 cm</< 10 14 cm. 

5.2. Magnetic field in W75N(B) 
5.2.1. Magnetic field strength 

The magnetic field strength along the line of sight was deter- 
mined from circular polarization measurements for 20 water 
masers. We detected significant (> 3cr) magnetic fields toward 
10 masers in VLA 1 and 7 masers in VLA 2. The detected fields 
along the line of sight (By) range from -400 mG to +600 mG for 
the H 2 masers around VLA 1 and from -200 mG to 1000 mG 
for those associated with VLA 2 (see Table [TJ. The wide range 
obtained for B^^ 1 and B^ 42 is not anomalous, and a similar 
range was reported for Cepheus A (Vlemmings et al. 2006). The 
changing in the sign indicates the reversal of the magnetic field 
(negative towards the observer, positive away from the observer) 
as already observed in other sources (e.g., Vlemmings et al. 
2006). Since the 6 values are close to 90°, i.e. close to the plane 
of the sky, a slight difference in these angles can produce an in- 
version of the sign of the magnetic fields. In the case of VLA 1, 
the positive and negative magnetic field strengths measured for 
VLA IB and VLA 1C groups indicate a twisted magnetic field, 
which is not the case for VLA 2. Apart from the H2O masers 
VLA2.01 and VLA2.02, which are located on the right edge of 
the ellipse 2 and show positive values, the other masers in the up- 
per left part of the ellipse 2 (from VLA2.09 to VLA2.26) present 
a clear separation between positive and negative magnetic field 
strengths (see Fig.|4]and Table[T]i. 

Based on these detections the absolute weighted magnetic 
field strengths along the line of sight, where the weights are 
Wi = 1/e 2 and e\ is the error of the i-th measurement, are 
<|fi^ LA1 |> = (81 ± 62) mG and (\B^ LA2 \) = (145 + 110) mG. In 
the case of VLA 1, it is more correct to determine for only 
the unsaturated masers; i.e., QBY^ 1 1) = (74 + 50) mG. The large 
errors are due to the large scatter of the magnetic field strength 
values, which depend on the projection effects, on the density of 
each water maser region, and on the local velocities of shocks. 
Because of this scatter it is difficult to evaluate the true strength 
of the full magnetic fields of the entire region. To compare the 
results obtained using different maser species we can estimate 
them roughly. Since <0) V lai = 83°^° D and <0> VL A2 = 85°+f° 
the absolute weighted magnetic field strengths are |Z?vlaiI = 
<|7i[ LA1 |>/«™<(9>vLAi ~ 700mG and |B V la2| ~ 1700mG, re- 
spectively. Such high values are also found by Vlemmings et al. 
(2006) in Cepheus A. As their shock velocity is 10 kms -1 , mea- 
suring a magnetic field of 600 mG, these authors explain this 
high value with the presence of a nearby magnetic dynamo. 

If the magnetic field is important throughout the collapse 
of a spherical cloud, the cloud forms primarily by flows along 
field lines, and the conservation of magnetic flux and of the 
mass imply B oc « 047 (Crutcher 1999). Although the scaling of 
the magnetic field in water maser is due to shocks, Vlemmings 
(2006, 2008 ) empirically showed that even water masers with 
number densities up to 10 n cirT 3 follow this relation, even if in 
non-masing gas of similar densities the magnetic field strengths 
are likely lower due to ambipolar diffusion. Hence, from the 
magnetic field strengths, we can determine the number den- 
sity of the cloud where the H2O masers arise by the relation 
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B cc n . Considering the pre-shock magnetic field (B\\ — 
16 mG) and density («h 2 = 10 9 cm 3 ) obtained by S09 from 
methanol maser emissions, we have «^ LA1 = 3 x 10 10 cirT 3 and 
w vla2 - i x io 1 1 cm 4 . Similar calculation can be done consid- 
ering the magnetic field strength obtained during the OH flare 
near VLA 2 (Slysh et al. 2010). They measured a magnetic field 
strength of about 70 mG, which is comparable to the total mag- 
netic field strength inferred by S09 (B = 50 mG). Since from 
the OH maser observations is generally possible to get only 
B and not By and since the flare occurred near VLA 2, in this 
case we have to use the value of |Bvla2| instead of the value 
of (|B|f LA2 |). For an OH number density of 10 8 cm 3 , we find 
n H^ A2 ~ 9 x 10 ll) cirT 3 for the water masers around VLA 2, close 
to what is obtained considering the methanol masers. These re- 
sults, which are close to the extreme upper limit of 10 11 cm -3 
of the water maser thermalisation (Elitzur et al. 1989), confirm 
the assumptions on the CH3OH and OH number densities and 
reinforce the measurements of the magnetic field strength made 
by S09, Slysh et al. (120101) and in this paper. The three indepen- 
dent magnetic field strength values obtained from three different 
maser species located at different positions show that the mag- 
netic field in massive star-forming regions can reach strength 
larger than expected. This large unexpected values of B might 
be due either to the presence of a magnetic dynamo or to the 
presence of very high-velocity shocks that strongly compress the 
gas, and consequently increase the magnetic field strength from 
50-70 mG (pre-shock methanol-OH region) to > 1 Gauss (wa- 
ter maser region). For the last hypothesis we refer the reader to 
§E3] 

With regard to the orientation of the magnetic field, it is 
worthwhile to discuss the influence of the Faraday rotation in 
our observations. 

5.2.2. Faraday rotation 

The measurements of the linear polarization angle (x) might be 
disturbed by the foreground, ambient, and internal Faraday rota- 
tion, which are given by 

<D[°] = 1.35 x 10~ 15 D [cm] n e [cm 4 ] B„ [mG] v~ 2 [GHz], (5) 

where D is the length of the path over which the Faraday ro- 
tation occurs, n e and By are respectively the average electron 
density and the magnetic field along this path and v the fre- 
quency. For the foreground Faraday rotation (i.e. the rotation due 
to the medium between the source and the observer) 0f is 0°4 at 
22 GHz, assuming the interstellar electron density and the mag- 
netic field are n e su 0.012 cm 4 and B\\ as 2 jjG, respectively. 

All water masers are not in the same plane, but they arise 
in clouds at different depths. This depth is impossible to deter- 
mine from our observations, so in order to estimate the ambient 
Faraday rotation (i.e. the rotation due to the ambient medium 
where masers arise), we assume the maximum projected separa- 
tion as upper limit of D among masers close to each other with 
the highest velocity difference, which is Z) max « 2 x 10 15 cm. 
For the electron density, we assume the value reported by Fish 
& Reid (2006), « e ~ 300crrT 3 , which can produce a rotation 
of 90° along the path amplification of OH maser around H 11 re- 
gions. Considering By = 16 mG (S09), the ambient Faraday ro- 
tation is <p. d < 26°. 

The internal Faraday rotation, which can destroy the linear 
polarization (Fish & Reid 120061 ). depends on the type of shock 
that pumps the water masers. They can either be a dissociative 



jump shock (J-shock) or a nondissociative continuous shock (C- 
shock). In the case of a C-shock, the intrinsic Faraday rotation at 
22 GHz can be considered negligible because the ionization state 
of the gas is controlled by cosmic -ray ionization, which gener- 
ates electron-molecular ion pairs at a rate of 10~ I7 s~' (Kaufman 
& Neufeld 1996 ), while we have to roughly estimate its contribu- 
tion for a J-shock. In this case, the electron density is n e = Xe w h 2 > 
where Xe is the ionization fraction (^ e = 10~ 5 - 10~ 4 ; Kylafis & 
Norman 1987 ) and «h, the number density of the H2O masers. 
In the most conservative situation, for which n e ~ 10 4 cm~ 3 , 
we determine an internal Faraday rotation of tens of degrees, 
enough to destroy the linear polarization in some cases. This 
means that either the electron density must be much lower or, 
more likely, the internal Faraday rotation must be negligible, and 
consequently the shock must be a C-shock. 

5.2.3. Magnetic field orientation 

Suppose that the water masers in both sources are pumped by C- 
shocks (see § 15.31 1. this means that the Faraday rotation is < 26°. 
Since the H2O masers associated with VLA 1 and VLA 2 indi- 
cate that there are two different magnetic fields around the two 
sources, it is worthwhile discussing them separately. All water 
masers show 9 > 8 a -i t ~ 55°; i.e., the magnetic field direction is 
perpendicular to the linear polarization vectors and close to the 
plane of the sky. 

The group VLA 1A shows an orientation of the magnetic 
field i£g LA1A = 70° +9°, which is in good agreement with what 
is reported by S09 (tp B = 73°+10°). This also confirms that the 
magnetic field orientation in this part of the source is close to 
the direction of the large-scale molecular bipolar outflow (66°) 
(S09). The other two groups VLA IB and C show disordered lin- 
ear polarization vectors, which cannot give information on the 
orientation of the magnetic field. The different orientation of the 
linear polarization vectors for these two groups might come from 
the different depths of the masers, i.e. high values of the ambient 
Faraday rotation, and therefore we can assume that the masers 
with angles far from -20° are located at higher depths than those 
with angles close to -20° (e.g. VLA1.14). 

Although around VLA 2 some H2O masers show the 90°- 
flip of the linear polarization angle (e.g., VLA2.15), which was 
also observed in some water masers in Cepheus A (Vlemmings 
et al. 2006), the magnetic field appears to be radial (see Fig.©. 
Considering (xi) ~ -12°, the magnetic field orientation is about 
18° that indicates a misalignment with the large-scale molecu- 
lar bipolar outflow and the magnetic field of VLA 1 , even though 
the inclination of both ellipses (~ 60°) is close to them. 

The different orientation of the magnetic fields of the two 
sources again implies different local physical environments for 
VLA 1 and VLA 2, and suggests that VLA 1 is the powering 
source of the large-scale molecular bipolar outflow, as already 
indicated by S09. 

5.3. J- or C-shocks? 

As reported in previous papers, the H2O masers are thought to 
be pumped by the transition of shocks generated in the two pro- 
tostars (e.g., Torrelles et al. 1997). These shocks can be either 
C-shocks or J-shocks. In a water masing region, C-shocks show 
v < 50 kms 4 and a post-shock T < 4000 K (Kaufman & 
Neufeld 1996), while J-shocks have velocities v > 45 kms -1 
and a post-shock plateau temperature that depends on the pre- 
shock density, the shock velocity, and the observed linewidth of 
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the water masers (Elitzur et al. 119891 ). Kylafis & Norman ( 1991) 
investigated the effects of the temperature on their J-shock model 
up to a temperature of about 1000 K, and they found that the 
brightness temperature starts to increase more slowly when a 
temperature -500 K is reached. 

We first consider VLA 1 . Assuming the proper motion of the 
water masers reported by T03, i.e. 19 kms -1 , as the velocity of 
the shock and as pre-shock density, the density of the methanol 
masers «h, = 10 9 crrT 3 (S09), the equation reported in Elitzur et 
al. ([19891 "for a J-shock 

rplateau * 40 °li0W* lOOkms-i lkms-ij K ' (6) 

where here Avvlai = 0.5 kms -1 , gives us a most likely post- 
shock temperature r^ t A1 u « 900 K, which is higher than the 
temperature measured considering the weighted intrinsic veloc- 
ity obtained from our fit (TVlai ~ 400 K). Moreover, from the 
model of a C-shock described by Kaufman & Neufeld ( 1996), 
we are able to estimate the velocity of the shock by considering 
the temperature 7Vlai ■ The model for a temperature of 400 K 
gives v^ LA1 ~ 15 km s _1 , which is close to the proper motion of 
the H2O masers. Therefore, following our arguments, the shock 
in VLA 1 is most likely a C-shock. 

In the case of VLA 2, the plateau temperature due to a J- 
shock is a 1000 K (where Av VL A2 = 0.7 kms" 1 , and 

vvla2 = 46 kms~'), which is less than the temperature that we 
measured ( 7Vla2 ~ 2500 K). The difference of the two temper- 
atures may be due to the turbulence velocity of the gas, which 
can be determined by considering AV\. In fact, AV; is related to 
the thermal and turbulence velocities by the equation 

W = V Ay t h + AV t 2 urb , (7) 

and considering the Eq.[3]we are able to estimate the turbulence 
velocity of the gas if a J-shock is present. The thermal velocity, 
which stems from the warming up of the gas by the compression 
due to the passage of a shock, for a temperature of 1000 K is 
AV t ^ LA2 a 1.6 kms -1 , so we have AV turb ~ 1.9 kms -1 . If in- 
stead we suppose that the shock is a C-shock from the model of 
Kaufman & Neufeld (fT996t for T WLA2 ~ 2500 K, we get v^ LA2 ~ 
43 kms _I , which is close to the expansion velocity that we de- 
termined in the present paper (vvla2 = 46 kms~'). Both veloc- 
ities are also in good agreement with the velocity of 43 km s -1 
measured by Slysh et al. ( 120101 ) for the OH masers near VLA 2 
obtained by multiepoch observations. However, vvla2 is close 
to the limit of 45 kms~', which is considered the limit between 
the two types of shocks, so we cannot rule one of them out by 
only considering the velocity. Since we obtain well-ordered high 
linear polarization fraction and linear polarization vectors, we 
have to suppose that the Faraday rotation is at most ~20°; oth- 
erwise, the vectors no longer appear to be aligned with the large 
structure (i.e., ellipse 2). This condition is only met if the in- 
ternal Faraday rotation is small; in other words, if the shock is 
a C-shock (Kaufman & Neufeld [1996 ). Therefore, we strongly 
suggest that there is a C-shock also in VLA 2 even if a J-shock 
cannot be completely ruled out. 
From the equation 

fi masei- _ „„ I "H 2 V's ,gs 

[mG] ~ V[!0 8cm ~ 3 ] [10km s" 1 ]' 

where «h 2 is the pre-shock density and v s the velocity of the 
shock, we can verify our assumption about the pre-shock den- 
sity and the rough estimation of |BvlaiI and |Bvla2|- Equation 



[8]is valid for both types of shocks (Kaufman & Neufeld Tl 9961 ). 
Considering «h 2 = 10 9 cirT 3 as before, we get |BvlaiI ~ 500 mG 
and |Bvla2| ~ 1200 mG, which are consistent with the rough val- 
ues reported in § 15.21 making the presence of a magnetic dynamo 
unnecessary. Finally, from our results we can strongly depict the 
following scenario for VLA 1 : the hot gas in a C-shock pumps 
the water masers and at the same time the warm dust associ- 
ated with this shock emits infrared photons, which pumps the 
methanol masers in the pre-shock region, as already suggested 
in S09. 



6. Conclusions 

We observed the 22 GHz water masers in full polarization mode 
in the massive star-forming region W75N with the VLBA. We 
detected 124 water masers around the two radio sources VLA 1 
and VLA 2, which appear to be in two different evolutionary 
stages; i.e., VLA 1 is more evolved and a radio jet has not 
formed yet in VLA 2. From our observations, we have shown 
that these two sources separated in the sky by ~ 1" (2000 
AU) have different local physical environments. In particular, 
the linear polarization indicates that the magnetic field around 
VLA 1 is aligned well with the large-scale molecular bipolar 
outflow, while the magnetic field is well-ordered around the 
shell-like expanding gas associated with VLA 2. The detection 
of Zeeman-splitting indicates an absolute weighted magnetic 
field strengths |B V laiI ~ 700 mG and |B V la2| ~ 1700 mG. 
We were also able to determine the type of shocks that pump 
the water masers associated with both sources. They both are 
C-shocks, even though a J-shock cannot completely be ruled out 
in the case of VLA 2. We suggest VLA 1 as the driving source 
of the large-scale molecular bipolar outflow. 
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